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Using a sample of 1.8 million DD mesons collected at the ?/>(3770) with the CLEO-c detector, we 
study the semileptonic decays D° — > -K~e + v e , D + — > 7r°e + ^ e , D° — > K~e + u e , and D + — » K°e + v e . 
For the total branching fractions we find B{D° ir~e + v e ) = 0.299(11)(9)%, B{D + -> -K°e + u e ) = 
0.373(22)(13)%, B{D° -» K~e+v e ) = 3.56(3)(9)%, and B(D+ -> K°e+v e ) = 8.53(13)(23)%, where 
the first error is statistical and the second systematic. In addition, form factors are studied through 
fits to the partial branching fractions obtained in five q 2 ranges. By combining our results with recent 
unquenched lattice calculations, we obtain |V cd | = 0.217(9) (4) (23) and |V CS | = 1.015(10)(11)(106), 
where the final error is theoretical. 



PACS numbers: 12.15.Hh, 13.20.Fc, 14.40.Lb 

I. INTRODUCTION 

In the standard model of particle physics, mixing of the 
quark mass eigenstates in their charged current interac- 
tions is described by the Cabibbo Kobayashi Maskawa 
(CKM) matrix Q. This 3x3 quark mixing matrix must 
be unitary and can be described by four independent pa- 
rameters. If the standard model is complete, experimen- 
tal determination of the CKM matrix elements should 
verify its unitarity. Deviations from unitarity would indi- 
cate the presence of physics beyond the standard model. 
A variety of CP-conserving and CP- violating observ- 
ables probe the elements of the CKM matrix and allow 
us to over-constrain it. Many of the key observables re- 
quire great precision or great sensitivity to provide the 
constraints at the level needed to test the validity of the 
standard model description. It thus remains a continuing 
experimental challenge to test the unitarity of the CKM 
matrix fully. 

Study of the semileptonic decay of D mesons plays 
a primary role in our understanding of the CKM ma- 
trix. These decays allow robust determination of the 
CKM matrix elements \V CS \ and \V c d\ by combining mea- 
sured branching fractions with form factor calculations, 
such as those based on unquenched lattice QCD (LQCD) 
0| ■ In addition, these measurements will provide preci- 



sion tests of LQCD itself [3j. One approach to tests of 
LQCD assumes unitarity of the CKM matrix and com- 
pares the constrained matrix elements Q to elements 
obtained with a combination of CLEO-c measurements 
and lattice form factors. A second approach, which is 
independent of CKM elements and thus free from the 
unitarity assumption, compares the measured and calcu- 
lated ratios of semileptonic and purely leptonic branching 
fractions. Verification of lattice calculations at the few 
percent level will provide validation for use of the lattice 
in the B system, where they are relied upon for several 
crucial theoretical quantities. 

This article presents a study of the D° — > n~e + v e , 
D+ -► TT°e+v e , D° -► K~e + v e and D+ K°e+v e de- 
cay modes (charge conjugate modes implied). A sum- 
mary of the analysis is also provided in a shorter com- 
panion article The results are based on a sample of 
1.8 million DD pairs collected with the CLEO-c detec- 
tor at the Cornell Electron Storage Ring (CESR) from 
281 pb" 1 of e+e" data at the -0(3770) resonance. The 
sample is a superset of, and approximately five times 
larger than, the data used to obtain the first CLEO-c 
semileptonic branching fraction measurements .6]. For 
each mode we determine the partial branching fractions 
in five q 2 ranges, with the sum of the five rates deter- 
mining the total branching fraction. Fits to the rates de- 
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termine the form factor shapes. By incorporating LQCD 
calculations into the form factor fits, we extract values 
for the CKM elements \ V c d\ and \ V CS \. Previous quenched 
lattice predictions carried errors of about 20%. Current 
unquenched LQCD calculations allow theoretical evalua- 
tion of the form factors at the 10% [2| level, with future 
improvement to the few percent level expected. 

Within this article, Section UT1 provides an overview of 
the formalism for exclusive semileptonic decays of charm 
mesons and their associated form factors. Sections Mil 
through I VII cover the experimental procedures for event 
reconstruction and extraction of the branching fractions, 
the systematic uncertainty evaluation, and the branching 
fraction results. Sections I VIII and I Villi explore the form 
factor shape constraints from our data and the extraction 
of | V cs | and \ V c d\- Section llXl presents our conclusions and 
comparisons with previous measurements. 



II. EXCLUSIVE CHARMED SEMILEPTONIC 
DECAYS 

The matrix element describing the semileptonic decay 
of a D meson to a pseudoscalar meson P is of the form 



M(D -> Pe+v e ) 



-i^V:„L»H u 



(1) 



where Gp is the Fermi constant, V cq is the appropriate 
CKM matrix element and and arc the lcptonic and 
hadronic currents. The leptonic current can be written 
in terms of the electron and neutrino Dirac spinors, u e 
and t>„, 



(2) 



In the case of pseudoscalar decays, where there is no 
axial-vector contribution, the hadronic current is given 
by 



H„ = (P(p)\q^c\D(p')} 



(3) 



where p' and p are the four-momenta of the parent D 
meson and the daughter P meson, respectively. The 
hadronic current is fundamentally a non-perturbative 
quantity that is difficult to evaluate. We can, however, 
re-parameterize the current by expressing it in terms of 
the independent four-momenta in the process, which for 
a pseudoscalar-to-pseudoscalar decay are the two four- 
momenta p' + p and q = p' — p. We can identify q as 
the four-momentum of the virtual W boson. A typi- 
cal formulation of the hadronic current in terms of these 
four-momenta is given by 



(p(p)\rfc\D(j/)) 



(4) 
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where Mp> is the mass of the D meson and mp is the 
mass of the final state pseudoscalar meson. The non- 
perturbative contributions are incorporated in the scalar 
functions f+(q 2 ) and foil 2 ), the form factors of the decay. 
Kinematic constraints require /+(0) = /o(0). A further 
simplification arises due to the small mass of the electron 
because q^L^ — > in the limit m e — > 0. Thus including 
only the /+ form factor in the hadronic current, 

{P{p)\qrc\D(p')) = f+(q 2 )(p'+PV, (5) 

is a very good approximation. With this form for the 
hadronic current the partial decay width becomes 



dT(D -> Pev e ) 
dq^ 



G F \V cq \ 
24tt 3 



2\|2 
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(6) 



The partial decay width (Eq. [5]) clearly reveals that 
extraction of the CKM matrix elements from measured 
rates requires prediction of the semileptonic form factors. 
Theoretical calculation of the form factors therefore has 
become a considerable industry. We focus here on pa- 
rameterizations of the form factors that we employ in our 
form factor studies and in extraction of \V c d\ and |V CS |. 

The goal of any particular parameterization of the 
semileptonic form factors is to provide an accurate, and 
physically meaningful, expression of the strong dynamics 
in the decays. To that end, one may express the form fac- 
tors in terms of a dispersion relation, an approach that 
has been well established in the literature (see for exam- 
ple Ref. Q and references therein). It is common to write 
the dispersive representation in terms of an explicit pole 
and a sum of effective poles: 



U(Q 2 ) 



/+(0) 



1 



Pk 



1 — a l — 



^ I i g* 

T fc wip olc 



(7) 



where pk and 7^ are expansion parameters. Given the 
underlying c — > q quark transition of the semileptonic 
decay, the mass m po i e is the mass of the lowest-lying cq 
vector meson. The parameter a gives the contribution 
from the vector meson pole at q 2 = 0. Using this disper- 
sion relation the true form factor can be approximated to 
any desired degree of accuracy by keeping sufficient terms 
in the expansion. This approach has the drawback that 
the decay dynamics are not explicitly predicted. Addi- 
tionally, experimental data have suggested the need for 
only a few parameters in the description of the form fac- 
tor shape. It is therefore natural to seek simplifications 
of this parameterization that can still capture the correct 
dynamics. 

Removing the sum over effective poles entirely, leaving 
only the explicit vector meson pole, provides one simplifi- 
cation route that is typically referred to as "nearest pole 
dominance" or "vector-meson dominance" . The result- 
ing "simple pole" parameterization of the form factor is 
given by 
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Experimental data disagree with the physical basis for 
this approximation, since measurements of the parame- 
ter TOp i c that fit the data do not agree with the expected 



vector meson masses 14] . Effectively, at low or medium 
values of q 2 the spectrum is distorted compared to the 
simple pole model, receiving contributions from the con- 
tinuum of effective poles above the lowest lying pole mass. 

The modified pole, or Becirevic-Kaidalov (BK) param- 
eterization Q, was proposed to address this problem. 
The parameterization keeps the first term from the ef- 
fective pole expansion, while making simplifications such 
that the form factor can be expressed using only two pa- 
rameters: the intercept /+(0) and an additional shape 
parameter. 1 The parameterization is typically expressed 
in the form 



Mo 2 ) 



/+(o) 



(i 



-)(i 



-) 



(9) 



This parameterization has recently been widely used in 
the extraction of semileptonic form factors from experi- 
mental measurements [3, E3, EM El . In addition, some 
recent LQCD calculations of the form factor have re- 
lied on this parameterization for extrapolation and inter- 
polation purposes 0, ■ This scheme requires several 
assumptions to reduce the multiple parameters initially 
present (Eq. [7]) to one. The BK ansatz assumes that the 
gluon hard-scattering contributions (5) are close to zero 
and that scaling violations (/3) are close to unity, which 
may be succinctly expressed as 



1 + 1/(3- 5 = 



(Ml 



/+(0) 



dq 2 



(10) 



q 2 =0 



Once again, however, the experimental data do not bear 
out these assumptions [3]. We should observe a ~ 1.75 
in order to obtain 1 + 1//3 — 5 = 2, whereas the observed 
data are removed from such values by many standard 
deviations. 

We note that both functional forms can provide ade- 
quate parameterizations of the data if their parameters 
are allowed to be non-physical. Without a physical un- 
derpinning for the parameterization, however, parame- 
ters obtained from theory and/or from different exper- 
iments may not agree if their form factor sensitivities 
differ as a function of q 2 . 

Our primary form factor shape analysis therefore uti- 
lizes a series expansion around q 2 = to that has been 
advocated by several groups for a physical description of 
heavy meson form factors d, [13, EL G3 • The series ex- 
pansion is congruous with the dispersion relations, and is 
guaranteed to contain the true form factor, yet is still rich 



1 There will be three parameters if the /o(<? 2 ) form factor, which 
we are neglecting due to the small electron mass, is also taken 
into account. 



enough to describe all variations that affect the physical 
observables. 

To achieve a convergent series, the expansion is formu- 
lated as an analytic continuation of the form factor into 
the complex t = q 2 plane. There is a branch cut on the 
real axis for t > {Mo + Mk,tv) 2 that is mapped onto the 
unit circle by the variable z, defined as 



z(q 2 ,to) 



\A+ - g 2 - yt + - tp 

y/t+ -<? + V*+ - *o ' 



(11) 



where t± = (Md ± tik,k) 2 and to is the (arbitrary) q 2 
value that maps to z = 0. The expression for the form 
factor becomes 

1 oo 

^^ = ^w^)S Gfc(io)W92,to)]fe ' (12) 



with 



P(q 2 ) 



1, D -> 7T 

z(q 2 ,M 2 D »), D^K 



(13) 



The P{q 2 ) factor accommodates sub-threshold reso- 
nances, which overcomes the convergence issues that a 
naive expansion would face with a nearby pole. Good 
convergence properties are expected since the physical 
region is restricted to \z\ < 1. The physical observables 
do not depend on the choice of 0(q 2 ,t o ), which can be 
any analytic function, or on the value of to- We report 
ak parameters that correspond to tg = and the "stan- 
dard" choice for (j> (see, e.g. Ref. pjj and Appendix [A}, 
which results from bounding ^2 °I from unitarity con- 
siderations. Appendix [Al presents results for an alternate 
choice of to that minimizes the maximum value of | z | over 
the physical range. If the series converges quickly, as ex- 
pected, it is likely that only the first two or three terms 
will be able to be seen in the data. We will explore the 
number of terms needed to adequately describe our data. 

While our primary form factor and CKM results will be 
based on the series expansion, we will also provide results 
based on the two pole parameterizations for comparative 
purposes. 



III. EVENT RECONSTRUCTION AND 
SELECTION 

The analysis technique rests upon association of the 
missing energy and momentum in an event with the 
neutrino four- momentum [20| , an approach enabled by 
the excellent hermeticity and resolution of the CLEO- 
c detector [3, [2l|. Charged particles are detected over 
93% of the solid angle by two wire tracking chambers 
within a 1.0 T solenoid magnet. The momentum reso- 
lution is 0.6% at 800 MeV/c. Specific ionization mea- 
surements from the tracking system in combination with 
a ring imaging Cerenkov detector (RICH) [22| provide 
particle identification. A CsI(Tl) crystal electromagnetic 
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calorimeter provides coverage over about 93% of 4ir, and 
achieves a typical ir° mass resolution of 6 MeV/c 2 . 

Electron candidates are identified above 200 MeV/c 
over 90% of the solid angle by combining information 
from specific ionization with calorimetric, RICH and 
tracking measurements. The identification efficiency, 
which has been determined from data, is greater than 
96% above 500 MeV/c and greater than 90% above 
300 MeV/c. The average probability that a hadron is 
misidentified as an electron is less than 0.8%. Below 300 
MeV/c the efficiency falls rapidly, reaching 60% in the 
200 — 250 MeV/c region. To reduce our sensitivity to 
final state radiation (FSR), we add photons within 3.5° 
of the initial electron momentum back into the tracking- 
based four-momentum. 

Charged pions and kaons from the signal decay are 
identified using specific ionization and RICH measure- 
ments. Pion candidates below 750 MeV/c and kaon 
candidates below 500 MeV/c are identified using only 
specific ionization information, which is required to be 
within three standard deviations (a) of that expected for 
the assigned particle type. For pion candidates above 
650 MeV/c, we also require the pion mass hypothesis be 
more likely than the kaon mass hypothesis. Above these 
momenta, candidate tracks must also pass RICH identi- 
fication criteria. Specifically, we require that pion (kaon) 
candidates are more than 3tr closer to a pion (kaon) hy- 
pothesis than a kaon (pion) hypothesis. 

A 7T° candidate must have a 77 mass within 2.5cr of the 
7T° mass. Kg candidates are reconstructed using a con- 
strained vertex fit to candidate n + n~ daughter tracks. 
The ir + w~ mass must be within 4.5a of the K s mass. 

To reconstruct the undetected neutrino we utilize the 
hcrmeticity of the CLEO-c detector to find the miss- 
ing energy and momentum. In the process e + e~ — > 
?A(3770) — ► DD, the total energy of the beams is im- 
parted to the DD system. Because the beam energies are 
symmetric and the beam crossing angle is small at CESR, 
each produced D has an energy of the beam energy to 
within a small correction. The missing four-momentum 
in an event is given by p miss = (E missi p miss ) = p tota] - 
DPchargcd - SPncutrai, where the event four- momentum 
Ptotai is known from the energy and crossing angle of the 
CESR beams. Charged and neutral particles for the sums 
must pass selection criteria designed to achieve the best 
possible |pmiss| resolution by balancing the efficiency for 
detecting true particles against the rejection of false ones. 

For the charged four- momentum sum, J^Pcharged, opti- 
mal selection is achieved with topological criteria. These 
criteria minimize multiple-counting that can result from 
low-momentum tracks that curl in the magnetic field, 
charged particles that decay in flight or interact within 
the detector, and spurious tracks. Tracks that are actu- 
ally segments of a single low transverse momentum "curl- 
ing" particle arc identified by selecting reconstructed 
track pairs with opposite curvature whose innermost and 
outermost diametric radii each match within 14 cm and 
whose separation in is within 180° ± 20°. For physics 



use we select the track segment that will best repre- 
sent the original charged particle based on track quality 
and distance-of-closest-approach to the beam spot. We 
employ similar algorithms to identify particles that curl 
more than once, creating three or more track segments. 
We also identify tracks that have scattered or decayed in 
the drift chamber, causing the original track to end and 
one or more tracks to begin in a new direction. We keep 
only the track segment with the majority of its hits be- 
fore the interaction point. Spurious tracks are identified 
by their low hit density and/or low number of overall hits 
and rejected. 

Each hadronic track must be assigned a mass hypoth- 
esis to calculate its contribution to the total energy sum. 
We assign a most probable mass hypothesis by combining 
detector measurement with particle production informa- 
tion. The production information is introduced because 
it is only statistically advantageous to identify a track 
as a kaon at a momentum where many more pions than 
kaons are produced when the detector's particle ID in- 
formation strongly favors a kaon. For each track, we 
first calculate a likelihood for the kaon and pion hypothe- 
sis based on specific ionization and RICH measurements. 
Those likelihoods are then weighted by the Monte Carlo 
(MC) prediction for the relative K~ and ir~ abundances 
in D decays at that track's momentum, which then gives 
us the true probability for each mass hypothesis. 

For the neutral four-momentum sum, ^Pneutrai, clus- 
ters resulting from the interactions of charged hadrons 
must be avoided. As a first step, calorimeter showers 
passing the standard CLEO proximity-matching (within 
15 cm of a charged track) are eliminated. Optimizations 
also revealed that all showers under 50 MeV should be 
eliminated. The processes that result in reconstructed 
showers ( "splitoffs" ) separate from but within about 25° 
of a proximity-matched shower tend to result in an en- 
ergy distribution over the 3x3 central array of the splitoff 
shower that "points back" to the core hadronic shower. 
We combine this information with the ratio of energies in 
the 3 x 3 to 5 x 5 arrays of crystals, whether the shower 
forms a good ir°, and the MC predictions for relative 
spectra for true photons versus splitoff showers to pro- 
vide an optimal suppression of the contribution. 

Association of the missing four-momentum with the 
neutrino four-momentum is only accurate if the event 
contains no more than one neutrino and if all true par- 
ticles are detected. For events with additional missing 
particles or doubly-counted particles, the signal modes 
tend not to reconstruct properly while background pro- 
cesses tend to smear into our sensitive regions. Hence, 
it is worthwhile to reject events for which independent 
measures indicate these problems. We therefore ex- 
clude events that have either more than one electron 
or non-zero net charge. Multiple electrons indicate an 
increased likelihood for multiple neutrinos, while non- 
zero net charge indicates at least one missed or doubly- 
counted charged particle. 

After application of the above criteria approximately 
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90% of the signal MC |p m i ss | distribution is contained in 
a central core with a ~ 15 MeV/c. 

To further enhance the association of the missing mo- 
mentum with an undetected neutrino in our final event 
sample, we require that the M^ iss = S^ iss — |p m iss| 2 be 
consistent with a massless neutrino. The M^ iss resolu- 
tion, 

cr(M 2 liss ) = 2E miss a(E miss ) © 2|p miss |cr(|p miss |), 

is dominated by the -E m i ss term since the resolution of 
IPmissI is roughly half that of -E m i ss - MC simulation in- 
dicated an optimal requirement of |M^ liss /2|p miss | | < 0.2 
GeV/c 3 , which (noting E miss w |p miss | for signal) pro- 
vides selection at approximately constant -E m i ss reso- 
lution. Additionally, because of the superior |p m iss| 
resolution, in subsequent calculations we take p v = 

( l^miss 1 5 Pmiss ) ■ 

Semileptonic decays D — ► Peu, where P is a pion or 
kaon, are identified by their consistency with the ex- 
pected D energy and momentum. Candidates are se- 
lected based on AE = (Ep+E e +E v ) — -Bbeam (expected 
to be zero within our resolution of about 20 MeV) and 
yields are extracted from the resulting distributions in 
beam-constrained mass Mb c (equivalent to D momen- 
tum and expected to be close to the known D mass). 
These quantities are corrected for the small boost re- 
sulting from the 3 mrad beam crossing angle. Because 
the \p v \ resolution dominates the AE resolution, we can 
improve our p v measurement by scaling it by the fac- 
tor ( satisfying (E P + E e + (E u ) - E hcam = 0. We 
use C,p v for the neutrino momentum in computation of 
M bc = \J El cam - \p P +p e + (pu\ 2 . The resulting reso- 
lution for M bc is 4 MeV/c 2 . 

Selection criteria were optimized by studying MC sam- 
ples independent of those used elsewhere in the anal- 
ysis. Sources of backgrounds include events with fake 
electrons, non-charm continuum production (e + e~ — ► qq, 
e + e~ — ► t + t~ and e + e~ — > jtp(2S)), and DD processes 
other than signal. 

The optimal AE requirement was determined to be 
-0.06 < AE < 0.10 GeV. For the Cabibbo-favored 
modes, the background level remaining after this se- 
lection is only a few percent of the signal level. For 
the Cabibbo-suppressed modes, there remains significant 
background from cross-feed among the signal modes, par- 
ticularly from the kaon modes, as well as from the re- 
lated modes D + — > K { [e + iy e and D + — > Kge + v e where 
Kg 7r + 7r~. Since the cross-feed typically involves par- 
ticles from the "other D" decay, we obtain some suppres- 
sion of this background with a g 2 -dependent requirement 
on A£ n s . for the non-signal particles in the event. We ob- 
tain A_E n s . by summing the energy of all non-signal parti- 
cles in an event, even though we do not specifically recon- 
struct the non-signal D decay. This criterion effectively 
imposes an additional constraint on the quality of the 
reconstructed neutrino. We also require D + — > 7r°e + ^ e 
candidates to have the smallest \AE\ compared to any 
other final state candidates in the event, and that these 



events contain no reconstructed D° — > K~e + v e candi- 
date. These criteria suppress cross-feed from the charged 
pion and kaon modes with almost no loss of true ir Q e + v e 
decays. The average background level (g 2 -dependent) in 
the pion modes is about 20% of the signal level. 

To simplify the statistical interpretation of our re- 
sults, as well as to suppress cross-feed from the Cabibbo- 
favored into the Cabibbo-suppressed modes, we limit the 
number of multiple entries per event such that a given 
event can contribute to at most one D° or one D + final 
state. For events with multiple D + candidates or mul- 
tiple D° candidates satisfying Mb c > 1.794 GeV/c 2 , we 
choose the candidate with the smallest |AS|, indepen- 
dent of q 2 . 

From the measured electron and the re-scaled neu- 
trino four-momenta we calculate q 2 = M 2 ^. from q 2 = 
[p v +p e ) 2 . The resulting resolution is 0.01 GeV 2 /c 4 , in- 
dependent of q 2 . 



IV. EXTRACTION OF BRANCHING 
FRACTIONS 

A. Method and Binning 

For each of the four signal modes we construct the 
Mbc distributions in five q 2 ranges: q 2 < 0.4 GeV 2 /c 4 , 
0.4 < q 2 < 0.8 GeV 2 /c 4 , 0.8 < q 2 < 1.2 GeV 2 /c 4 , 
1.2 < q 2 < 1.6 GeV 2 /c 4 and q 2 > 1.6 GeV 2 /c 4 . These 20 
distributions are fit simultaneously to extract the partial 
branching fraction for each interval. The total branch- 
ing fraction is then obtained in each mode by summing 
its five partial branching fractions. Fitting in five q 2 
ranges minimizes the experimental sensitivity of the to- 
tal branching fractions to form factor shape uncertainties, 
while simultaneous fitting of all four modes ensures self- 
consistent handling of the cross-feed backgrounds among 
the modes. 

The fit utilizes a binned maximum likelihood approach 
extended to include the finite statistics of the MC sam- 
ples following the method of Barlow and Beeston [26| . 
The Mbc distribution is divided into fourteen uniform 
bins over the range 1.794 < M bc < 1-878 GeV/c 2 . 



B. Fit Components and Parameters 

We fit the data to the signal components and five back- 
ground components. The signal mode components are 
obtained from MC generated using EvtGen [2J] and mod- 
ified pole-model form factors l8j with parameters from 
the most recent LQCD results 121 . W e apply several cor- 
rections to our GEANT-based [23| MC samples to im- 
prove simulation of the neutrino reconstruction proce- 
dure. 

From independent studies, mostly based on CLEO- 
c samples with one fully reconstructed hadronic D de- 
cay, we evaluate corrections and associated systematic 
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uncertainties for simulation of hadronic showers, false 
charged particles and charged-particle identification. We 
find that the simulations of charged particles, charged- 
particle momentum resolution and photon-energy reso- 
lution need no correction, though we include the uncer- 
tainties in the systematic uncertainty evaluation. We 
reweight the MC samples to correct the rate and spec- 
trum for Kl production (which affects the neutrino- 
reconstruction efficiency), for ir° and ir~ production 
in our full D decay model, and for the momentum- 
dependent rate at which a K~ fakes a 7r~. All of these 
corrections affect the cross-feed background rates into 
and between the Cabibbo-suppressed modes. They lead 
to few percent (or less) changes in the measured yields, 
but are determined to better than 10% of themselves. 

In the MC samples we select only true electrons (recon- 
structed tracks that have been matched to a generator- 
level electron) with a probability for acceptance given by 
data-measured efficiencies described earlier. We thereby 
exclude from the MC any events caused by identification 
of a fake electron and instead estimate this background 
using data, as we describe in detail below. This proce- 
dure eliminates any reliance on MC predictions for cither 
electron efficiency or the rate at which hadrons fake elec- 
trons. 

We are sensitive to the distortion of efficiency and 
kinematics in our signal modes due to FSR. Based on 
the angular and energy distributions for FSR photons, 
we correct our signal MC, generated with the PHO- 
TOS package without interference effects included, 
to the Kaon Leading-Order Radiation (KLOR) [28| cal- 
culations modified for charm decay. 

For each reconstructed q 2 interval in a given mode, we 
generate a MC sample in the same (generator level) q 2 
interval, to which the full analysis is applied. That is, we 
obtain the full set of 20 reconstructed Mb c distributions 
from each of these 20 independent samples. For each of 
the generated q 2 intervals, a single floating parameter, 
which corresponds to the efficiency-corrected data yield 
in that interval, controls the normalization of all its 20 
reconstructed distributions. The relative normalizations 
among those reconstructed distributions remains fixed at 
the level predicted by our corrected MC. Because the sig- 
nal rate in each reconstructed range drives the normal- 
ization for the corresponding generated q 2 interval, the 
data in effect fixes the cross-feed rates into the other 19 
reconstructed distributions. 

We also use MC samples to describe the DD back- 
ground and the three continuum contributions. We ab- 
solutely scale the continuum components according to 
their cross sections at the ^(3770) and the measured 
data luminosity. The non-signal DD sample was gen- 
erated using EvtGen, with decay parameters updated 
to reflect our best knowledge of D meson decays. This 
component floats separately for each reconstructed final 
state, but the relative rates over the five q 2 regions within 
that state are fixed. This approach helps to reduce our 
sensitivity to inaccuracies in the D decay model. Fi- 



nally, we input MC components for D + — » K\e + v e and 
D + — > Kg(Tr°Tr°)e + i' e , whose rates in each q 2 region are 
tied to those for the signal D + — ► K < g(jT + Tr~)e + v e mode. 

The contributions from events in which hadrons have 
faked the signal electron are evaluated using data. The 
momentum-dependent electron identification fake rates 
from pions and kaons are measured using a variety of data 
samples. We obtain our background estimates by ana- 
lyzing a data sample with no identified electrons. Each 
track in each event in this sample is treated in turn as the 
signal electron. The contribution in each mode is then 
weighted according to the fake rate. The fake electron 
component is then added to the fit with a fixed, abso- 
lute, normalization. 

Finally, we allow the fit to adjust the M^ c resolution in 
the D° -> n-e + v e , D+ -> 7r°e+^ e , and L>° -> K~e + v e 
modes by applying a Gaussian smear to these distribu- 
tions. As a result the signal MC Mb c resolution in these 
modes is increased from ^3.5 MeV/c 2 to match the data 
resolution of ~4 MeV/c 2 . The M^ c resolution in the 
D + — > K*ge + v e signal MC matches the data resolution 
very well so we apply no additional smearing to this mode 
in the fit. 

In summary we have 27 free parameters in the fit: the 
20 signal rates, the 4 non-signal DD normalizations and 
the 3 Mbc smearing parameters. This leaves us with a 
total of 280 — 27 = 253 degrees of freedom for the fit. 



C. Checks and Results 

The direct fit results are displayed as plots in Mb c , di- 
vided into the appropriate q 2 ranges. These are shown 
in Fig. [TJ The value of the likelihood for this fit is 
— 21n£ = 275.5 for 253 degrees of freedom. Note that 
each of the 20 distributions is described by mainly one 
free parameter - the signal rate within that bin, with a 
more constrained contribution from the parameters (reso- 
lution and DD background rate) that float independently 
for each signal mode, with the relative contribution into 
each q 2 interval fixed within a mode. All other contri- 
butions are either explicitly or effectively fixed by other 
constraints. 

Other important reconstructed kinematic variables are 
presented integrated over q 2 with the components scaled 
according to the nominal fit. Fig. [5^ shows the AE dis- 
tributions for events within the signal-enhanced region 
\M hc -M D \ < 0.015 GeV/c 2 . Fig.^b shows cos6 We , the 
cosine of the angle between the W in the D rest frame 
and the electron in the W rest frame, in the signal Mb c 
and AE regions. All of our signal modes should exhibit 
a sin 2 Owe dependence independent of the form factor 
though acceptance effects distort the reconstructed dis- 
tribution. The fits describe the observed distributions 
very well. Finally we find that our fit generally agrees 
well with the observed momentum (p e ) spectrum for the 
signal electron (Fig. [2t). The poorest agreement is ex- 
hibited by the Tr°ev e mode, where the probability of x 2 
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FIG. 1: M bc distributions for the modes (a) D° -> %~e + u e , (b) D + -» 7r°e+^ e , (c) D° -> K~e + v e , and (d) D + -y K%e + v e . 
The data are shown by the points, and the fit components (histograms) are normalized using the nominal fit results (see text): 
signal MC (clear), cross- feed and non-signal DD MC (gray), continuum MC (light gray) and e + fakes (black). 
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TABLE I: Results of fit to DD MC sample with statistics 
of 40 x £data for all q 2 bins. Yinput is the true yield, Vat the 
efficiency corrected yield from the fit, and (TY rit the la error 
on the efficiency-corrected fit yield. 







( 


Yinput — Vflt)/0Y fit 










True 


q 2 interval (GeV 2 /c 4 ; 


i 




Decay 


< 0.4 


0.4-0.8 


0.8 - 1.2 


1.2 - 1.6 > 1.6 


All q 2 




0.55 


-0.92 


0.98 


-0.33 


1.16 


0.51 


7v°e+u e 


-1.57 


0.37 


-0.55 


0.85 


0.98 


-0.95 


K~e + v e 


-2.34 


1.27 


1.54 


0.18 


0.99 


-0.14 


K° s e+u e 


-0.29 


-0.49 


1.77 


1.14 


0.54 


0.76 



is still over 3%. 

To test the fitting procedure, we fit a set of mock data 
with known input branching fractions created from the 
large DD MC sample (~40 x£data) used to obtain our 
non-signal DD background estimate. We fit the sam- 
ple using distributions from our standard signal MC and 
from the non-signal portion of the generic DD sample. 
Because our "data" in this case derives from the same 
underlying decay model and detector simulation as our 
fit inputs, we do not apply the corrections noted in the 
previous section that remove data/MC differences. 

Table U presents the differences between our measured 
and the generator- level rates. We see no biases at this 
greater level of sensitivity, demonstrating the reliability 
of the fitting procedure. Furthermore, the semilcptonic 
modes in the generic MC sample used to simulate the 
"data" were generated using ISGW2 [2^] form factors, 
which have a significantly different q 2 behavior than the 
LQCD-derived form factors 2] of our signal MC. Our test 
therefore also verifies that we have adequately subdivided 
the q 2 range to avoid significant dependence on input 
form factor modeling, even at levels significantly more 
sensitive than we can probe with the current data. 



V. EXPERIMENTAL SYSTEMATIC 
UNCERTAINTIES 

The systematic uncertainties for the D° modes are 
summarized in Tables [TT1 and IIII1 The first table presents 
the complete list, while the second breaks down the 
neutrino-reconstruction simulation errors into compo- 
nent parts. The corresponding systematics tables for the 
D + modes are presented in Tables IIVI and [V] For indi- 
vidual uncertainties we give the sign of the error relative 
to the change in the lowest q 2 range. The largest system- 
atic uncertainties are those associated with the number 
of DD pairs (needed for normalization in the branching 
fraction determination, as described in Section VI) and 
with neutrino reconstruction simulation. Uncertainties 
in neutrino simulation include both inaccuracies in de- 
tector simulation and uncertainty in the decay model of 
the non-signal D, as discussed above. 

The starting point for the assessment of many of the 



systematic uncertainties is the measurement of any dis- 
crepancies between data and MC in the desired quantities 
(e.g., signal pion efficiency, signal kaon efficiency, etc.). 
Such measurements (or limits) are made using an inde- 
pendent data sample - in most cases, a sample of events 
with one of the two D mesons from the VK3770) fully 
reconstructed in a hadronic mode. In the case of signifi- 
cant discrepancies, the MC samples are corrected for use 
in our nominal fit (the fit used to obtain our final branch- 
ing fraction results, as opposed to any of the fits used to 
obtain systematic uncertainties) as noted above. Such 
corrections lead to changes in the measured yields of up 
to a few percent, but are determined precisely enough 
to yield sub-percent systematic uncertainties. For each 
systematic category, we determine the size of its contribu- 
tion by biasing the MC samples away from their nominal 
configuration at the level given by the uncertainty of the 
independent study. We re-fit the data with these biased 
MC samples, and use the deviation of the fit results from 
their nominal values to provide an estimate of the uncer- 
tainty. We note that because of the correlations among 
the five q 2 intervals in a given mode, the sum over q 2 
of the systematic errors tends to be less sensitive to the 
systematic variations than the individual intervals them- 
selves. 

The number of DD pairs, used to convert the mea- 
sured yields to branching fractions (see below), is a di- 
rect product of the CLEO-c hadronic branching fraction 
analysis [3(3] . We combine the statistical and systematic 
uncertainties from that analysis for our uncertainty esti- 
mates. 

We have assessed the uncertainties associated with the 
finding and identification efficiency for each of the signal 
hadrons. For the signal and tt^, the charged track- 
finding efficiency is already accounted in the tracking ef- 
ficiency portion of the v simulation uncertainty. They 
have, however, additional particle identification (PID) 
criteria associated with them, for which we assess a cor- 
rection and uncertainty. For the signal ir° and Kg we 
assess a correction and uncertainty for the reconstruc- 
tion efficiencies of these particles. We evaluate each of 
these four uncertainties by first measuring a momentum 
dependent, and hence q 2 dependent, correction and fit 
the measurements with a linear parameterization. The 
best fit result is applied as a correction in the nominal fit. 
To evaluate the systematic uncertainty, we identify the 
largest systematic variation on the x 2 = 1 ellipse from 
the linear fit. The branching fractions are most affected 
by the largest variation in overall normalization, while 
the form factors (Section IVII[) arc most affected by the 
largest variation in slope. 

We have uncertainties associated with the electron 
identification efficiency and the rates for hadrons to mis- 
reconstruct as (fake) electrons. We vary the efficiency 
and fake rates used in the analysis of our MC samples 
(see Section HVBp according to the uncertainties from 
the data studies used to measure them, and re-fit the 
data to evaluate our sensitivity. 
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FIG. 2: The kinematic distributions for (a) AE, (b) cosdwe, and (c) p e , for events falling within the M^ c signal region for 
each of the four signal modes. The data are shown by the points, and the fit components (histograms) are normalized using 
the nominal fit results (see text): signal MC (clear), cross-feed and non-signal DD MC (gray), continuum MC (light gray) and 
e + fakes (black). The dotted lines in (a) indicate for each mode the AE region used in fitting. 



Modeling of tt* and ir° production - spectra and rates 
- in D decay significantly affects the background shape 
and rate for the signal cross-feed background into the 
pion signal modes. The large effect results because a 
pion from the non-signal D decay can be swapped in as 
the signal pion candidate. We measure the background 
pion spectra in data using the inclusive D decays on 
the "other side" of a fully reconstructed hadronically- 
decayed D "tag", and correct the MC spectra accord- 
ingly. To be conservative in the associated systematic 



uncertainty, we take the full difference for results ob- 
tained using corrected and uncorrected spectra. In the 
signal D° — > ir~e + i/ e mode we also correct the cross-feed 
background from D° — *■ K ~e + v e events that results from 
misidentifying a as a rr^ . Once again the uncertainty 
estimate is taken as the difference of our measured rates 
obtained using the corrected and uncorrected fake rates. 

For three of our signal modes, D° — ► it~e + v e , D° — > 
K~e + v e and D + — > n°e + v e , we have systematic uncer- 
tainty associated with the additional Mb c resolution pa- 
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TABLE II: Summary of full and partial branching fraction systematic errors (%) for the D° — ► ir~e + v e and D° — ► A'^e+^e 
signal decay modes. The sign represents the direction of change relative to the change in the [0, 0.4) GeV 2 /c 4 interval in each 
mode. 



q 2 interval (GeV 2 /c 4 ) q 2 interval (GeV 2 /c 4 ) 



Systematic 


[0,0.4) 


[0.4,0.8) [0. 


8, 


1.2) 


[1.2, 


1.6) 


[l-6,9max] 


All 


<1 2 


[0,0.4) | 


0.4,0.8) 


[0.8, 


1.2) 


[1.2, 


1.6) 


[1-6, g max ] 


All 


<f 


Number DD 


1.51 


1.51 


1 


.51 


1 


.51 


1.51 


1 


.51 


1.51 


1.51 


1 


.51 


1 


.51 


1.51 


1 


.51 


u e simulation 


1.45 


1.77 


2 


.21 


2 


.87 


1.59 


1 


.69 


1.52 


1.99 


1 


.96 


2 


.39 


1.28 


1 


.80 


tv° efficiency 


0.01 


0.01 





.01 





.00 


-0.03 


-0 


.00 


0.00 


0.00 





.00 





.00 


0.01 





.00 


Kg efficiency 


0.01 


-0.01 


■0 


.02 


-0 


.04 


-0.08 


-0 


.03 


0.01 


0.00 





.00 





.00 


0.07 





.01 


7T" PID 


0.44 


0.41 





.39 





.39 


0.40 


0. 


41 


0.00 


0.00 





.00 





.00 


0.00 





.00 


K~ PID 


0.05 


0.02 





.02 





.01 


0.02 





.03 


0.32 


0.30 





.28 





.26 


0.26 





.30 


e+ PID 


0.76 


0.40 





.75 





.36 


0.44 


0. 


.56 


0.65 


0.58 





.60 





.52 


0.46 





.61 


e + Fakes 


2.50 


0.45 


•0 


.01 





.05 


0.92 





,88 


0.57 


0.07 


-0 


.04 





.01 


0.39 





.25 


7r° production 


0.01 


0.02 





.03 





.01 


0.06 





03 


0.00 


0.00 


-0 


.01 





.00 


0.01 





.00 


7T _ production 


0.07 


0.42 





.44 





.10 


-1.96 


-0. 


.24 


0.00 


0.00 





.00 





.05 


0.23 





.01 


K~ fakes 


0.67 


1.01 





.93 





.35 


-0.08 





.58 


0.00 


0.00 


-0 


.01 


-0 


.01 


-0.02 





.00 


■K~e + u e Mhc res. 


0.93 


1.06 





.85 





.83 


1.05 


0. 


.95 


0.02 


0.01 





.01 





.01 


0.06 





.01 


K~e + v e Mbc res. 


0.07 


0.04 





.02 





.01 


0.01 





.03 


0.09 


0.08 





.07 





.09 


0.16 





.09 


7r"e + ^ e Mbc res. 


0.01 


0.02 





.00 





.01 


-0.13 


-0 


02 


0.00 


0.00 





.00 





.00 


0.00 





.00 


e + veto 


0.05 


0.04 


•0 


.01 


-0 


.14 


-0.01 


0. 


.00 


0.08 


0.07 





.00 


-0 


.02 


-0.03 





.05 


FSR 


0.85 


1.53 





.97 





.91 


0.75 





99 


0.63 


0.61 





.56 





.48 


0.47 





.59 


Model dep. 


0.50 


-0.01 


•0 


.09 





.43 


-1.55 


-0 


19 


0.33 


-0.11 


-0 


.16 


-0 


.41 


-1.29 





.02 


Total 


3.70 


3.26 


3 


.26 


3 


.56 


3.74 


2 


.95 


2.44 


2.66 


2 


.63 


2 


.95 


2.51 


2 


.53 



TABLE III: Summary of full and partial branching fraction systematic errors (%) associated with neutrino modeling in the 
MC for the D ' — ► -K~e + v e and D° -> K~ e+ v e signal decay modes. The sign represents the direction of change relative to the 
change in the [0, 0.4) GeV 2 /c 4 interval in each mode. 









D° 












D° -» K~e+v e 












q 2 interval 


(GeV 2 /c 4 


) 








q 2 interval 


(GeV 2 /c 4 


) 




v Systematic 


[0,0.4) 


[0.4,0.8 


) [0.8,1.2) 


[1.2,1.6) 


[l-6,<7max 


] All q 2 


[0,0.4) 


[0.4,0.8 


) [0.8,1.2) 


[1.2,1.6) 


[l-6,<?max] 


All q 2 


split-off showers 


0.58 


0.90 


1.92 


2.58 


0.48 


1.17 


0.89 


1.54 


1.49 


2.00 


0.67 


1.29 


K% showers 


0.15 


0.13 


0.11 


0.06 


-0.60 


-0.04 


0.01 


-0.02 


-0.04 


-0.05 


-0.06 


-0.01 


Kl production 


0.67 


0.66 


0.66 


0.67 


0.65 


0.66 


0.66 


0.65 


0.66 


0.67 


0.69 


0.66 


track efficiency 


0.60 


0.56 


0.50 


0.63 


0.45 


0.54 


0.37 


0.39 


0.44 


0.44 


0.21 


0.39 


track resolution 


0.00 


1.02 


0.01 


0.46 


0.85 


0.46 


0.28 


0.43 


0.50 


0.60 


0.34 


0.39 


split-off rejection 


0.58 


-0.02 


0.16 


0.04 


-0.22 


0.12 


0.59 


0.56 


0.48 


0.43 


0.21 


0.54 


particle ID 


0.01 


-0.02 


0.08 


0.02 


0.09 


0.04 


0.06 


0.04 


0.05 


0.03 


0.04 


0.05 


shower resolution 


0.03 


0.09 


0.06 


-0.01 


0.11 


0.06 


0.00 


0.00 


0.00 


0.00 


-0.01 


0.00 


fake tracks 


0.76 


0.71 


0.70 


0.71 


0.72 


0.72 


0.72 


0.72 


0.71 


0.71 


0.71 


0.72 


Total 


1.45 


1.77 


2.21 


2.87 


1.59 


1.69 


1.52 


1.99 


1.96 


2.39 


1.28 


1.80 



rameter. The statistical uncertainty already has a con- 
tribution from allowing this parameter to float. We esti- 
mate the contribution to the systematic uncertainty for 
each mode by increasing the value of that mode's reso- 
lution parameter by one standard deviation beyond the 
best fit result. 

We must also account for any uncertainty associated 
with modeling event loss from the single electron veto 
because of secondary electrons from photon conversions 
and other processes. According to data studies (using 
the CLEO-c "tagged" samples, where one of the two D 
mesons from the ■0(3770) is fully reconstructed), our MC 
simulation models the number of secondary electrons in 
our events accurately within the error of the study. The 



most likely potential source of uncertainty arises from 
mis-modeling the rate for photon conversion within the 
detector material. For the uncertainty estimate we there- 
fore vary this contribution over the range allowed by the 
maximum allowed uncertainty of our data study, about 
8%. 

For the systematic uncertainty associated with the fi- 
nal state radiation (FSR) modeling, we take the dif- 
ference between the KLOR and PHOTOS predictions. 
This simulates a change in the radiative branching frac- 
tion of up to 16% in the most extreme case. Because 
the majority of the correction results from the lack of 
the FSR interference terms between the charged hadron 
and electron, the systematic should be an overestimate of 
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TABLE IV: Summary of full and partial branching fraction systematic errors (%) for the D + — > ir°e + i/ e and D + — » _R'°e + i/ e 
signal decay modes. The sign represents the direction of change relative to the change in the [0, 0.4) GeV 2 /c 4 interval in each 
mode. 



D+ -> ir°e + u e D+ -» if°e + ^ e 

<3 2 interval (GeV 2 /c 4 ) <7 2 interval (GeV 2 /c 4 ) 



Systematic 


[0,0.4) 


[0.4,0.8) [0. 


8, 


1.2) 


[1.2, 


1.6) 


[1-6, <? ma x] 


All 


<f 


[0,0.4) 


[0.4, 


0.8) 


[0.8, 


1.2) 


[1.2, 


1.6) 


[1-6, <? max 


] All 


<f 


Number DD 


1.60 


1.60 


1 


.60 


1 


.60 


1.60 


1. 


.60 


1.60 


1 


.60 


1 


.60 


1 


.60 


1.60 


1 


.60 


v simulation 


2.54 


3.41 


2 


.57 


2 


.53 


2.45 


1. 


.96 


1.71 


1 


.75 


1 


.82 


1 


.84 


2.18 


1 


.74 


7T° efficiency 


0.87 


0.56 





.77 


1. 


.07 


1.07 


0. 


.85 


0.00 


0. 


.01 


-0 


.01 


-0 


.02 


0.00 





.00 


K% efficiency 


0.02 


0.02 





.08 





.11 


0.18 


0. 


.07 


1.05 


1. 


.00 





.94 





.88 


0.81 


1 


.00 


TV' PID 


0.17 


0.37 





.13 


-0 


.24 


-0.29 


0. 


.06 


0.01 





.01 





.00 


-0 


.03 


-0.02 





.00 


PID 


0.17 


0.37 





.12 


-0 


.24 


-0.29 


0. 


.06 


0.01 


0. 


.01 





.00 


-0 


.02 


0.00 





.00 


e+ PID 


1.13 


0.56 





.33 





.98 


0.01 


0. 


.62 


0.62 


0. 


.65 





.52 





.59 


0.76 





.61 


e + Fakes 


1.52 


0.14 


-0 


.29 


-0 


.07 


0.64 


0. 


.44 


0.38 


-0. 


.03 


-0 


.17 





.09 


1.00 





.14 


7T° production 


0.43 


0.81 





.76 


-0 


.73 


-1.87 


-0. 


.04 


0.02 


-0. 


.01 





.00 


-0 


.11 


-0.14 


-0 


.01 


7r~ production 


0.07 


-0.02 





.03 





.02 


1.46 


0. 


.29 


0.02 


-0. 


.01 





.00 





.19 


1.12 





.04 


if - fakes 


0.01 


-0.01 





.04 





.07 


0.20 


0. 


.05 


0.01 


0. 


.01 


-0 


.03 


-0 


.07 


-0.16 


-0 


.01 


7r _ e + i/ e Afbc res. 


0.00 


0.00 - 


-0 


.01 





.05 


-0.29 


-0. 


.05 


0.01 





.01 





.00 





.04 


0.09 





.01 


7^ _ e + i/ e M^c res. 


0.00 


0.01 





.00 





.00 


0.01 


0. 


.00 


0.02 





.01 





.00 





.00 


0.02 





.01 


7r°e + ^ e Mbc res. 


2.62 


1.27 


3 


.77 


1 


.17 


1.97 


2. 


.09 


0.01 


0. 


.01 





.01 





.05 


0.08 





.01 


e + veto 


0.26 


-0.01 





.20 





.03 


-0.14 


0. 


.07 


0.02 


0. 


.09 





.02 





.12 


-0.30 





.04 


FSR 


0.26 


0.48 





.47 





.68 


0.65 


0. 


.49 


0.25 


0. 


.46 





.55 





.64 


0.60 





.41 


Model dep. 


0.56 


0.08 


-0 


.08 





.76 


0.08 


0. 


.28 


0.35 


-0. 


.16 


-0 


.28 


-0 


.83 


-1.51 


-0 


.06 


Total 


4.57 


4.19 


5 


.00 


3 


.76 


4.52 


3. 


.53 


2.70 


2. 


.70 


2 


.73 


2. 


.87 


3.69 


2 


.67 



TABLE V: Summary of full and partial branching fraction systematic errors (%) associated with neutrino modeling in the 
MC for the D + — > Tv°e + i/ e and D + — » K°e + v e signal decay modes. The sign represents the direction of change relative to the 
change in the [0, 0.4) GeV 2 /c 4 interval in each mode. 









D+ - 


ir°e+v e 












D+ 














q 2 interval 


(GeV 2 /c 4 


) 










q 2 interval 


(GeV 2 /c 4 ) 






v Systematic 


[0,0.4) 


[0.4,0.8 


) [0.8,1.2) 


[1.2,1.6) 


[1-6, <?max] 


All 


? 2 


[0,0.4) 


[0.4,0.8 


) [0.8,1.2) 


[1.2,1.6) [1 


■6? 9max] 


All q 2 


split-off showers 


0.62 


2.94 


2.11 


1.30 


-1.68 


1. 


.14 


0.17 


0.37 


0.21 


-0.13 


-1.36 


0.18 


K% showers 


0.19 


0.17 


0.08 


0.13 


-0.83 


-0. 


03 


0.01 


0.06 


-0.08 


-0.15 


-0.46 


-0.02 


Kl production 


1.10 


1.07 


1.07 


1.08 


1.13 


1. 


.09 


1.07 


1.07 


1.09 


1.10 


1.09 


1.08 


track efficiency 


0.51 


0.37 


0.18 


-0.14 


0.13 


0. 


24 


0.62 


0.57 


0.70 


0.66 


0.39 


0.62 


track resolution 


1.08 


0.05 


-0.09 


-0.90 


0.13 


0. 


.12 


0.43 


0.45 


0.64 


0.49 


0.92 


0.49 


split-off rejection 


1.66 


1.08 


0.68 


1.45 


-0.81 


0. 


.86 


0.81 


0.88 


0.84 


0.99 


0.21 


0.84 


particle ID 


0.09 


-0.02 


0.04 


0.07 


-0.03 


0. 


03 


0.01 


0.02 


0.00 


0.01 


0.08 


0.01 


shower resolution 


0.00 


-0.01 


0.18 


0.30 


0.04 


0. 


08 


0.00 


-0.02 


0.03 


-0.01 


-0.09 


0.00 


fake tracks 


0.77 


0.71 


0.70 


0.71 


0.71 


0. 


72 


0.72 


0.70 


0.69 


0.69 


0.70 


0.71 


Total 


2.54 


3.41 


2.57 


2.53 


2.45 


1. 


96 


1.71 


1.75 


1.82 


1.84 


2.18 


1.74 



the FSR uncertainty from final or initial-state particles. 
This overestimate compensates for the unknown direct 
(structure-dependent) contributions. 

The final systematic error we assess is the dependence 
on our modeling of the form factor input to our signal 
MC. We reweight each of our signal MC samples with a 
different form factor input, namely ISGW2. The nominal 
form factor input to our signal MC is a BK parameteriza- 
tion [8( with parameters determined by lattice QCD [2|. 
The q 2 spectra of the latter differ markedly from those 
of ISGW2. We fit with the re-weighted MC spectra and 
the difference to the nominal fit gives the systematic error 
associated with model dependence. The small uncertain- 
ties obtained in this study confirm our conclusion drawn 



from fitting the large MC sample. 



VI. BRANCHING FRACTION RESULTS 

Combining the results of the fit and the systematic un- 
certainty estimates gives us the final efficiency-corrected 
yield measurement for each mode. From that yield 
(y), we obtain the branching fraction B — Y/2N D ]j, 
where N D ^ is the number of neutral (N D o[)o) or charged 
{N D + D -) pairs in our sample. We obtain these num- 
bers from an independent CLEO-c analysis [3(| based 
on the comparison of events with one reconstructed D 
to events with both D decays reconstructed, in certain 
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TABLE VI: Summary of the efficiencies (e) and efficiency-corrected yields for each q 2 interval and the corresponding partial 
branching fractions, the total branching fractions, the branching ratios and the isospin ratios. In all cases the first errors 
are statistical and the second are systematic. For the K° mode, the efficiency and yields correspond to the reconstructed 
Kg — > 7r + 7r~ mode, so do not include the initial production amplitude or 7r + 7r~ branching fraction factors. 

q 2 interval (GeV 2 /c 4 ) 







< 0.4 


0.4-0.8 


0.8 - 1.2 


1.2 - 1.6 


> 1.6 


Total 


£(%) 

Yield 


>e)(%) 


19.4 

1452(113)(49) 
0.070(5)(3) 


21.0 

1208(102)(35) 
0.059(5)(2) 


Z? u — ► 7r-e + ^ e 
22.4 

1242(99) (36) 
0.060(5)(2) 


22.8 

906(85)(29) 
0.044(4) (2) 


22.4 

1357(103)(46) 
0.066(5)(2) 


0.299(H)(9) 


e(%) 
Yield 

B(7T° e + 


u e )(%) 


7.5 

1379(168) (59) 
0.084(10)(4) 


8.0 

1584(180)(61) 
0.097(H)(4) 


D+ -» 7r u e+j/ e 
7.9 

1012(154)(48) 
0.062(9)(3) 


7.2 

1028(158)(35) 
0.063(10)(2) 


5.7 

1101(174)(47) 
0.067(H)(3) 


0.373(22)(13) 


e (%) 
Yield 
B(K~e 


+ Ve){%) 


19.2 

29701(441) (569) 
1.441(21)(35) 


20.5 

21600(377) (473) 
1.048(18)(28) 


D u K~e + v e 
20.0 

14032(304) (301) 
0.681(15)(18) 


18.3 

7001(225)(178) 
0.340(11)(10) 


13.9 

991(112)(20) 
0.048(5)(12) 


3.557(33)(90) 


e(%) 
Yield 
B(K°e 


V Ve){%) 


11.7 

19480(466) (417) 
3.436(82)(93) 


12.3 

14422(415)(306) 
2.544(73) (69) 


D+ -> K u e+u e 
12.5 

9009(327) (194) 
1.589(58)(44) 


12.2 

4656(236)(107) 
0.821(42)(24) 


12.5 

789(104)(26) 
0.139(18)(5) 


8.53(13)(23) 


Ro(%) 
R+(%) 

In 

Ik 




4.89(39)(12) 
2.45(31)(9) 
2.12(31)(9) 
1.06(3)(3) 


5.59(48) (12) 
3.80(45) (13) 
1.54(22) (7) 
1.04(4) (3) 


8.85(74)(15) 
3.89(61)(17) 
2.47(43)(13) 
1.09(5)(3) 


12.9(13)(2) 
7.6(12)(2) 
1.78(32) (7) 
1.05(6) (4) 


137(19)(3) 
48(10)(2) 
2.48(45) (13) 
0.88(15)(3) 


8.41(32)(13) 
4.37(27)(12) 
2.03(14)(8) 
1.06(2)(3) 



hadronic modes. For the same data set that we have 
used, that analysis finds N D o Bo = (1.031 ± 0.016) x 10 6 
and N D+D - = (0.819 ± 0.013) x 10 6 . Our fit yields, ef- 
ficiencies and branching fractions for each mode, in each 



q range, are presented in Table rVTl The total branching 
fractions for each mode (also listed in Table IVl|) are 



B{D° -> ir-e + is e ) = 0.299(11)(9)%, 



(14) 



B{D+ -> 7r°e+i/ e ) = 0.373(22)(13)%, (15) 



B(D° -> K~e + v e ) = 3.56(3)(9)%, 



and 



B{D+ 



K°e+u e 



= 8.53(13)(23)%. 



(16) 



(17) 



The errors listed are statistical and systematic, respec- 
tively. 

We also measure the branching fraction and partial 
width ratios in each q 2 range. The full results are given in 
Table I VII To determine the partial width ratios we used 
the Particle Data Group lifetimes [4] rrjo = 410.3 ± 1.5 
fs and td+ = 1040 ± 7 fs. For the integrated q 2 ranges 
we find the ratios of branching fractions 



B(D° 



B{D° -► K-e+v e ) 



.41(32)(13)% (18) 



and 



B(D+ -> TT°e + v e ) 
B{D+ -> K°e+u e ) 



= 4.37(27) (12)%. (19) 



The partial width ratios, which are expected to satisfy 
isospin relationships, are found to be 



.03(14)(8) (20) 



and 



Ik = 



T(D° - 


-> 7T e + v e ) 


T(D+ - 


-> TT°e + V e ) 


T(D° - 


■4 K-e + u e ) 


T(D+ - 





1.06(2)(3). (21) 



We expect I v = 2 and Ik — 1, hence the measured par- 
tial width ratios satisfy isospin symmetry within our ex- 
perimental precision. 



VII. FORM FACTORS 

For each of our four signal decay modes we have ob- 
tained partial branching fraction results in five q 2 ranges. 
To extract information about the form factors we use the 
relationship 



_ G 2 F \V cq \ 2 
1 24n 3 T D 



(22) 
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TABLE VII: Summary of form factor results for the series parameterization and pole model fits. Correlation coefficients for the 
total uncertainty between variables in any two (three) preceding columns are given by p (pij)- The first errors are statistical 
and the second are systematic. The values for the mode are isospin corrected. For the series parameters (aj) we have 

assumed \V CS \ = 0.976 and |V cd | = 0.224. 

Series Parameterization - Three Parameter Fits 



Decay 


a 


ai 


«2 


P01 P02 


pl2 


(0) 


1 + 1/(3-5 


P 


X 2 /d-o.f 


n~e + v e 


0.044(2) (1) 


-0.18(7)(2) 


-0.03(35)(12) 


0.81 0.71 


0.96 


0.140(7) (3) 


1.30(37) (12) 


-0.85 


2.0/2 


n°e+u e 


0.044(3) (1) 


-0.23(H)(2) 


-0.60(57)(15) 


0.80 0.67 


0.95 


0.138(H)(4) 


1.58(60) (13) 


-0.86 


2.8/2 


K-e+Ve 


0.0234(3) (3) 


-0.009(21)(7) 


0.52(28)(6) 


0.62 0.56 


0.96 


0.747(9) (9) 


0.62(13)(4) 


-0.62 


0.2/2 


K°e+u e 


0.0224(4) (3) 


0.009(32)(7) 


0.76(42)(8) 


0.72 0.64 


0.96 


0.733(14)(11) 


0.51(20)(4) 


-0.72 


1.7/2 








Series Parameterization 


- Two Parameter Fits 








Decay 


ao 


ffli 


P 


\V cq \f+(0) 




l + l/(3-S 


P 


X 2 /d.o.f 




n~e + u e 


0.044(2) (1) 


-0.173(19)(7) 


0.66 


0.140(5)(3) 




1.27(H)(4) 


-0.80 


2.0/3 




7r°e + u e 


0.046(2) (1) 


-0.124(30)(9) 


0.69 


0.147(7) (4) 




1.01(16)(5) 


-0.78 


4.0/3 




K~e + u e 


0.0230(2) (3) 


-0.047(6) (3) 


0.34 


0.734(6) (9) 




0.86(4) (2) 


-0.43 


3.8/3 




K°e+u e 


0.0218(3)(3) 


-0.046(9)(4) 


0.53 


0.713(9)(11) 




0.87(6)(3) 


-0.60 


4.9/3 








Simple Pole Model Fits 






Modified Pole Model Fits 




Decay 


l^c,|/+(0) 


m polc (GeV/c 2 ) 


P 


X 2 /d-o.f 




l^c«|/+(0) 


a 


P 


xV d.o.f 


n~e + u e 


0.146(4)(2) 


1.87(3)(1) 


0.63 


3.11/3 




0.142(4)(2) 


0.37(8) (3) 


-0.75 


2.1/3 


n°e+u e 


0.149(6)(3) 


1.97(7) (2) 


0.65 


4.42/3 




0.147(7) (4) 


0.14(16)(4) 


-0.75 


4.07/3 


K~e + u e 


0.735(5) (9) 


1.97(3)(1) 


0.36 


2.67/3 




0.732(6) (9) 


0.21(5)(3) 


-0.42 


4.32/3 


K°e + u e 


0.710(8) (10) 


1.96(4)(2) 


0.53 


4.1/3 




0.708(9) (10) 


0.22(8)(3) 


-0.59 


5.3/3 



to relate the form factor f+(q 2 ) to the partial branching 
fraction B in a particular q 2 range. In this expression, 
Td is the total decay width of the parent D meson, and 
i denotes the particular q 2 interval. A specific functional 
form is chosen for f+(q 2 ) (see Section|TT| and the parame- 
ter values are determined via a \ 2 fit to the five measured 
Bi. In order to account for the correlations between the 
branching fractions in each q 2 range we minimize the ex- 
pression 

X 2 = E( B *-w)Ci7 1 ( B i-fi)> ( 23 ) 

ij 

where yi is the fit prediction for the branching fraction in 
the i th q 2 interval, and is the inverse of the covari- 
ance matrix. The integration in each bin is performed 
numerically on each fit iteration using the trapezoidal 
rule. 

The systematic uncertainties on the form factor pa- 
rameters are evaluated using the same method as for the 
branching fraction analysis. We take the set of branch- 
ing fractions that result from the branching fraction fit 
for each systematic uncertainty, then redo the fit for the 
form factors. The difference in these fit parameters from 
the nominal results is taken as the estimate of the sys- 
tematic uncertainty. The list of systematic uncertain- 
ties evaluated is the same as for the branching fraction 
analysis (see Section [V}. Note that for the systematic 
errors found by exploring a one standard deviation x 2 
ellipse in normalization versus g 2 -dependence and tak- 
ing the largest observed deviation from the nominal re- 
sult (tt~ PID, K~ PID, tt° finding, and K% finding), the 
form factors and the branching fractions will have their 
largest deviations in very different regions of the ellipse. 
The form factor will be most sensitive to the region of 



the ellipse that causes the largest variation as a function 
of q 2 , while the branching fraction is most sensitive to 
the overall normalization. 

Fitting with the full covariance matrix that includes 
both statistical and systematic uncertainties and corre- 
lations (see Appendix iBj) yields almost identical central 
values and total errors. 

We evaluate the form factor shape using the functional 
form given by the series parameterization as described in 
Section [TTJ For comparative purposes we also provide 
results based on the two pole models described in Sec- 
tion [Til For the series model we perform fits using both 
the first two and the first three expansion parameters a^. 
This tests both our sensitivity to the number of param- 
eters in the expansion and the convergence of the series. 
We express our results in terms of the physical observ- 
ables, the intercept |V^g|/+(0) and 1 + 1/(3 — 6, as well 
as giving the expansion parameters. In the simple pole 
model we fit for the intercept and the pole mass m po i e , 
while in the modified pole model we fit for the intercept 
and the shape parameter a, which summarizes the ef- 
fective pole contribution. The results for all modes are 
summarized in Table IVHI Comparisons of the four fits, 
for each of the four modes, are shown in Fig. [3] To allow 
systematic differences to be viewed clearly between the 
various parameterizations, we normalize the data and all 
fit results to the result of the three parameter series fit 
in each q 2 interval. 

For the series expansion, comparison of the two- 
parameter and three-parameter fits shows that our kaon 
data prefer a non-zero quadratic z term. The probabil- 
ity of x 2 improves from 29% (22%) to 89% (44%) go- 
ing from two to three terms in the series for the K~ 
(K°) fit. The pion measurements currently lack sensitiv- 
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3600307-016 

• Data Simple Pole Modified Pole Series{2) Series(3) 



£ i-s 
< 

< 

1 



D u ^7i e + v P 



D + -> 7i e + v P 



D°-> K"e + v P 



0.5 1 1.5 2 2.5 0.5 1 1.5 

q 2 (GeV2/c 4 ) 

FIG. 3: Form factor fit comparison for all modes. All data 
(points) and fits (histograms) are normalized to the relevant 
three-parameter series fit result (Series(3), line at 1). The 
simple pole, modified pole, and two-parameter series fit (Se- 
ries(2)) are shown by triple-dot-dash, dashed, and solid his- 
tograms, respectively. 



ity to probe this term, and two and three parameter fits 
yield similar results for the first two parameters. Since 
a quadratic term appears to be preferred for the kaons, 
however, we include that term in our series fits to the 
pion data to improve the probability that our shape un- 
certainties bracket the true form factor shape. While the 
central value for ai is an order of magnitude larger than 
the other terms, we stress that regions of parameter space 
with <22 of similar magnitude to ao and at fall will within 
the 90% hypercontour for the fit, so no strong statements 
can be made about the size of ai or about the conver- 
gence (or potential lack thereof) of the series from these 
data. 

For the pole models we observe that the parameteri- 
zations can provide a shape that describes our data ad- 
equately, but only with parameter values that do not 
support their physical basis. Although the fits give quite 
reasonable x 2 values (see Table IVlT|) . the poles masses do 
not agree with the ¥n* (Md* ) masses expected for the 
kaon (pion) modes by over 3a for the most precise fits. 
The 1 + 1/P — S results from the D° — > K~e + v e series 
expansion fit are over 3a from the value of ^2 necessary 
for physical validity of the BK parameterization, while 
those derived from our a values for the kaon modes are 
tens of a away. 



VIII. EXTRACTION OF |V„| AND |V cd | 

We extract \V c d\ and \ V CS \ by combining our |T4 g |/+(0) 
results from the three parameter series expansion fits 
with the unquenched LQCD results Q (0) = 

0.64(3) (6) and f ( + D ^ K) (0) = 0.73(3) (7). For the D° -> 



7T and Z? H 



7T° modes we find \V C 



0.005 ± 0.023 and \V cd \ = 0.216 ± 0.017 ± 0.006 ± 0.023, 
respectively. For the D Q — > K~ and D + — > K° modes, 
we find \V CS \ = 1.023 ±0.013 ±0.013 ±0.107 and \V CS \ = 
1.004 ± 0.020 ± 0.015 ± 0.105. Averaging the D° and D+ 
results (taking into account correlated and uncorrelated 
uncertainties) we find 



\V cd \ =0.217 ±0.009 ±0.004 ±0.023 



and 



\V CI 



1.015±0.010± 0.011 ±0.106. 



(24) 



(25) 



The uncertainties, statistical, systematic and theoretical, 
respectively, are dominated by the discretization uncer- 
tainty in the LQCD charm quark action, which should 
be improved in the near future for the Fermilab action, 
or greatly reduced through the use of other actions.. 

We can also extract the ratio | V c d | / 1 V cs | from 
the ratio of our measured form factors. From 
the z expansion fits to our D° data, we ob- 



v cd \fi D ~ 



tain 
0.010 ± 0.003 



0/ 



(0)/\V cs \fi D ^ K \0) = 0.187 ± 
while from our D data we obtain 

|Kd|/j^WlK s |/+ D " K) (0) = 0.188 ±0.015 ±0.004. 
The errors are statistical and systematic, respectively, 
and all correlations have been taken into account. Aver- 
aging, again with correlated uncertainties accounted for, 
we obtain 



\Vcd\f. 



(D- 

+ 



>7r) 



(0) 



\v cs \fi D ^ K \o) 



0.188 ± 0.008 ± 0.002. 



(26) 



We can combine this result with calculations of 
/(0)^ x V/(0) (jD ^ K) to obtain the ratio of CKM ele- 
ments. A recent light cone sum rules (LCSR) calcula- 



tion, for example, obtains 
0.84 ± 0.04, which implies 



If 



J + 



\Vcd\ 

IF,., I 



= 0.223 ± 0.010 s t a t ± 0.003 syst ± 0.011 LCS r. (27) 



rd 



0.218 ± 0.011 ± 



IX. SUMMARY 

In summary, we have measured branching fractions 
and branching-fraction ratios for four semileptonic D de- 
cay modes in five q 2 bins. The branching fraction results 
are the most precise measured to date and agree well 
with world averages Q. Our modified pole a parameter 
results agree within 1.3<7 with previous determinations 
by CLEO III 0, FOCUS [l^, and Kev results from 
Belle[l8j], but show over 3a disagreement with Belle K\xv 
results and LQCD fits. The a parameters obtained with 
our individual Kev results are separated from the recent 
BaBar result [lj| by about 2.5<r. The z expansion re- 
sults between BaBar and our Kev agree closer to the 2a 
level or better, depending on the level of correlation be- 
tween the BaBar r\ and T2 parameters. The discrepancy 
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with LQCD is difficult to quantify because the covari- 
ance matrix for the LQCD form factors is lost during the 
chiral extrapolation procedure for the published analysis 
0. We have made the most precise CKM determina- 
tions from D semileptonic decays to date, and the re- 
sults agree very well with neutrino based determinations 
of \V c d\ and charmed-tagged W decay measurements of 
\V CS I \M ■ Overall, these measurements represent a marked 
improvement in our knowledge of D semileptonic decay. 
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Council of Canada. 



APPENDIX A: THE ^-EXPANSION: DETAILED 
FORMS AND ALTERNATE RESULTS 

The standard choice for the outer function (f>(t, t ) in 
the z expansion for f+(q 2 ) fEq. [T2")) arises from considera- 
tions of unitarity. From a perturbative Operator Product 
Expansion (OPE) calculation, one can show [^, [T^, [H| 
that the choice 

<Kt,to) = <* (y/t+ - t + y/t+ - t Q ) (Al) 

t+-t (yrp^ + yt + -t-) 3 / 2 
x (t+-to)V4 {ju—t + ^Y 

leads to a constraint on the coefficients 

n a 

£«2<1, (A2) 

for any choice of n a ■ The bound corresponds to forbid- 
ding the production rate of Dn states by the relevant 
current to exceed the inclusive production rate, which 
can be calculated within the OPE. To leading order, the 
coefficient a is given by 

Q =V^. (A3) 

Numerically, we have taken the charm quark mass to be 
m c = 1.2 GeV. 



The choice of the parameter to within the ^-expansion 
provides a potential source of ambiguity when compar- 
ing experimental results. In our fits, we have, for sim- 
plicity, chosen to = in our form factor fits utiliz- 
ing the z-expansion. Another common choice for to is 
that which minimizes the maximal value of the map- 
ping z(q 2 ) over the entire physical range. The value 

to = t+ ^1 — \/T— t-/t+) , where t± = mu ± rriK^ ac- 
complishes this minimization. The best fit a% values for 
our three parameter fit using this alternate value for to 
are presented in Table IVlTTI The values for |V C9 |/ + (0) 
(q = s,d) that we find in these fits are identical, within 
the precision we are quoting, to those presented in Ta- 
ble ECO 

Finally, some experimental results for the z-expansion 
are presented in terms of the ratios u = ai/ao for i > 0. 
To allow straightforward comparison, we also quote our 
results in this form in Table IVlTTl 



APPENDIX B: CORRELATION MATRICES 

To allow complete external use of the partial branch- 
ing fractions presented in this paper, we present the sta- 
tistical and systematic uncertainty correlation matrices. 
These matrices will allow, for example, for simultaneous 
fits of these results with other experimental results to ob- 
tain form factor parameters. The statistical correlation 
matrix (Table IIX[) is derived from the 20 x 20 covariance 
matrix produced in our fitting procedure. 

To obtain the systematic correlation matrix (Table |X|) . 
we create a separate covariance matrix from the corre- 
lated motions of all 20 yields in each individual system- 
atic study. We then sum the resulting matrices to obtain 
the total systematic covariance matrix. In the absence 
of correlations, this procedure would reduce to adding 
the systematic contributions for a given measurement in 
quadrature. In producing the covariance matrix for the 
form factor systematic uncertainty, we assume that the 
two pion modes are fully correlated and similarly for the 
two kaon modes, but treat the pion and kaon uncertain- 
ties as uncorrelated. For the N D + D - and N DO £,o un- 
certainties, we take into account the 39% correlation in 
those yields. 
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TABLE IX: The statistical correlation matrix obtained from the simultaneous fit to the data (see Section [IV B [I . The lines indicate the mode boundaries. The modes 
are labeled by their final state hadron. Within each submode, the five q 2 intervals are ordered from lowest to highest. 
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TABLE X: The total systematic correlation matrix for the 20 measured mode / q 2 intervals (see Section [V}. The lines indicate the mode boundaries. The modes are 
labeled by their final state hadron. Within each mode, the five q 2 intervals are ordered from lowest to highest. 
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